Optical properties of correlated materials — 
Generalized Peierls approach and its application to VO2 
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The aim of the present paper is to present a versatile scheme for the computation of optical 
properties of solids, with particular emphasis on realistic many-body calculations for correlated 
materials. Geared at the use with localized basis sets, we extend the commonly known lattice 
"Peierls substitution" approach to the case of multi-atomic unit cells. We show in how far this 
generalization can be deployed as an approximation to the full Fermi velocity matrix elements that 
enter the continuum description of the response of a solid to incident light. We further devise an 
upfolding scheme to incorporate optical transitions, that involve high energy orbitals that had been 
downfolded in the underlying many-body calculation of the electronic structure. As an application 
of the scheme, we present results on a material of longstanding interest, vanadium dioxide, VO2. 
Using dynamical mean-field data of both, the metallic and the insulating phase, we calculate the 
corresponding optical conductivities, elucidate optical transitions and find good agreement with 
experimental results. 

PACS numbers: 



I. INTRODUCTION 

Correlated matter is characterized by an enormous sen- 
sitivity with respect to changes of external parameters. It 
is the merit of this responsiveness that a remarkable rich- 
ness of properties emerges in these systems. Correlation 
effects seem, for instance, to be a vital issue to outstand- 
ing phenomena such as high temperature superconduc- 
tivity and colossal magnetoresistance. In the latter case, 
the possibility of tuning the fundamental behavior of a 
material by an external field will undoubtedly lead to 
yet improved data storage devices. A better understand- 
ing of the various effects of strong correlations is thus a 
highly desirable goal of condensed matter physics, both 
from the theoretical and the technological point of view. 

On the experimental side, numerous techniques have 
been devised for and applied to the study of correlated 
materials of ever growing complexity. Optical spec- 
troscopy, which is the subject of this paper, is, in a way, 
the most natural among them : Optical detectors are 
sampling the response to incident light, as do our eyes, al- 
beit accessing frequencies, and thus phenomena, that are 
beyond our vision. The technique is particularly suited 
for tracking the evolution of the system under changes of, 
for instance, temperature or pressure. This is owing to 
a generally high precision, and the fact that, contrary to 
e.g. photocmission spectroscopy or x-ray experiments, re- 
sults arc obtained in absolute values. Especially, the exis- 
tence of sum-rules allows for a quantitative assessment of 
transfers of spectral weight upon changes of the system 
properties. Therewith optical spectroscopy is particu- 
larly adapted to the study of correlated materials^ 3 -^. 

On the theory side, while weakly correlated mate- 
rials are well described within density functional the- 
ory (DFT)£, e.g. in the local density approximation 



(LDA}£, and moderate correlation effects are cap- 
tured by perturbative approaches, such as Hedin's GW 
approximation^, it was the advent of dynamical mean- 
field theory (DMFT)(for reviews see e.g^^), and its 
realistic extension, LDA+DMFT (for reviews se o 11 ' 12 ), 
which allowed for the description and understanding of 
several metal-insulator transitions that are derived from 
the Mott-Hubbard or related mechanisms. Though our 
discussion of optical properties within realistic many- 
body approaches is quite general, and applicable to other 
techniques, we present results on VO2 that are based on 
LDA+DMFT calculations. 

The paper is organized as follows : After having now 
expanded on the potency of optical spectroscopy and the 
motivation for more theoretical efforts, we will in the re- 
mainder of Section I briefly review experimental and the- 
oretical knowledge about vanadium dioxide. In Section 
II we develop our formalism for the optical conductiv- 
ity within realistic calculations. Section III is devoted 
to a detailed discussion on Fermi velocities. This part 
contains our major innovations. Readers less interested 
in technical details are welcome to jump directly to Sec- 
tion IV, which presents our theoretical optics spectra for 
vanadium dioxide in both, the metallic and the insulating 
phase. 



A. Vanadium dioxide — The material 

1. Basic electronic structure 

At its metal- insulator transition^ 3 . (T c =340 K), V0 2 
transforms from a metallic high temperature phase of ru- 
tile structure into an insulating monoclinic (Ml) phase, 
in which the vanadium atoms pair up to form tilted 
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dimers along the c-axis. This Ml phase was found to 
be non-magnetic*". 

Over the decades several scenarios were evoked to ex- 
plain the metal-insulator transition. This is not the em- 
phasis of the current paper, and we shall only briefly 
summarize the basic electronic structure. 

VO2 has a vanadium 3d 1 configuration and the crystal 
field splits the 3d-manifold into t2 S and empty e^ com- 
ponents. The former are further split into e^ and ai g 
orbitals. In the rutile phase these orbitals overlap, ac- 
counting for the metallic character. In the Ml phase, the 
&i g split into bonding/anti-bonding orbitals, due to the 
aforementioned vanadium dimerization. 

The Goodenough scenario 15 ' 16 of the insulator advo- 
cates the structural effect of the unit-cell doubling due 
to the dimer formation as the main origin of the gap- 
formation, and thus attributes the insulating behav- 
ior to a Pcicrls transition^ 7 -. Zylbersztejn and Motfei 8 - 
on the other hand stressed the importance of local 
Coulomb interactions, and thought the transition to be 
of, what we call today, the Mott-Hubbard type. Exper- 
iment o 19 ' 20 ' 21 ' 22 ' 23 ' 24 were interpreted to support one or 
the other of the two scenarios. Important to note is that 
neither of the two phases are well-described within stan- 
dard band-theory approaches. In the rutile phase these 
miss bandwidth-narrowing and satellite features, as seen 
e.g. in photoemissior*2i, and the bad metal conductiv- 
ity seen in transport measurements^ 5 -. In the Ml phase, 
the problem is even more fundamental, since band- 
theory fails to produce an insulating behavio r 26 i 27 ' 28 . 
For reviews see e.gi 29 ' 30 . As we will detail later, the 
LDA+DMFT approach succeeds in describing the ex- 
perimental findings of both, the metalli c - 31 ' 32 ' 33 and the 
insulating^ 3 , phase. As a matter of fact the current 
calculations of optical properties rely on our previous 
LDA+DMFT work^i^ which in particular extended 
on the interpretation of the nature of the insulating Ml 
phase. Compatibility with experimental results on the 
optical conductivity strengthens this picture. For the in- 
sulating phase, also GW type of calculation o 30 ' 36 ' 37 ' 38 as 
well as LDA+U based approache a 28 ' 39 open the charge 
gap. Further, VO2 has also been studied within cluster 
based methods, see e.g . 40 ' 41 ' 42 . 



2. Insights by optical measurements 

Optical measurements on VO2 were first performed by 
Barker et aL— , and Verleur et a/.—. By probing different 
orientations of single crystal samples, they evidenced an 
anisotropy in the optical response of the Ml insulator. 
More precisely, the conductivity depends on whether the 
electric field is parallel or perpendicular to the crystal- 
lographic rutile c-axis when going below the transition 
temperature. This is to be expected from the changes 
in the crystal-structure and the unit-cell doubling along 
the c-axis. This anisotropy was confirmed by ultraviolet 
reflectance measurements^ and x-ray experiment s 21 ' 46 



(see also^ 7 -). Ladd et al^- performed experiments under 
pressure, and noticed that c-axis stress reduces the tran- 
sition temperature considerably more than is the case 
for hydrostatic pressure. Okazaki et alX studied re- 
flectance spectra of thin films with an orientation of the 
electric field perpendicular to the rutile c-axis as a func- 
tion of temperature, and found indications for electron- 
phonon coupling. Recent studies by Qazilbash et a/4 8 - 
(see also 49 ' 49 ) on polycrystallinc films with preferential 
[010] orientation 5 ^ confirmed the bad metal behavior of 
rutile VO2 evidenced in transport experiments^ 5 -. Indeed, 
rutile VO2 is found to violate^ 8 - the Ioffe-Regel-Mott limit 
for resistivity saturation 5 ^, i.e. the electron mean free 
path is comparable to, or smaller than the lattice spac- 
ing and Boltzmann transport theory breaks down. As 
we shall see below, pronounced differences in the optical 
response are found between the individual experiments. 

Theoretically, the optical response of Ml VO2 was in- 
vestigated by means of a self-consistent model GW cal- 
culation by Continenza et al£^, which was found to im- 
prove on LDA results for the dielectric function, when 
comparing with experiments^. Also, a clear polariza- 
tion dependence was evidenced. 

Further, the dielectric response of both the metallic 
and the insulating phase were calculated within LDA by 
Mossanek and Abbate^ 3 -. In the metallic phase, peak 
positions and the polarization dependence were qualita- 
tively captured. The issue of the bad metallic behav- 
ior was not addressed, which is natural since it lies be- 
yond band theory. As to the insulating Ml phase, a 
rigid shift was introduced to the LDA band-structure, 
such as to "artificially" produce a gap. This proce- 
dure, again, resulted in qualitative agreement with ex- 
periment. However, we believe that the electronic struc- 
ture is characterized by an enhanced &\ g bonding/anti- 
bonding splitting 3 ^, which is not reproduced by an 
orbital- independent shift. An orbital-dependent onc- 
particle potential, on the other hand, actually docs cor- 
rectly capture spectral properties to a surprising de- 
gre o 34 ' 35 . We will come back to this in Section III. 



II. OPTICAL CONDUCTIVITY IN REALISTIC 
CALCULATIONS 

A. Optical conductivity from DMFT calculations 

Within the field of strongly correlated electrons, cal- 
culations of the optical conductivity within the DMFT 
framework were first performed by Jarrell et al£*, and 
Pruschke et ai£^ for the case of the Hubbard model. 
Rozenberg et al^-^ studied the phenomenology of 
the different optical responses of the Hubbard model 
throughout its phase diagram in conjunction with exper- 
iments on V2O3. In the realistic LDA+DMFT context, 
optical conductivity calculations were first performed 
by Bliime r 58 ' 59 for the case of degenerate orbitals. A 
more general approach, was developed in£°- for the study 
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of transport properties. Further, recent LDA+DMFT 
works that use simplified approaches to the Fermi veloc- 
ities can be found in^^i. Our work goes along the lines 
of the mentioned approaches. We will however use a full 
Hamiltonian formulation, therewith allowing for the gen- 
eral case of non-degenerate orbitals, and we extend the 
intervening Fermi velocities to multi-atomic unit cells, 
which becomes crucial in calculations for realistic com- 
pounds. 

Alternative techniques were presented by Perlov et 
al£2- in the Korringa-Kohn-Rostoker (KKR) context, and 
by Oudovenko et al*&. The idea in the latter work is to 
diagonalize the interacting system, which allows for the 
analytical performing of some occurring integrals due to 
the "non-interacting" form of the Green's function. Ow- 
ing to the frequency-dependence of the self-energy, how- 
ever, the diagonalization has to be performed for each 
momentum and frequency separately, so the procedure 
may become numerically expensive. 

First accounts of the presented optics scheme have 
been given in Refp^ for V2O3 (based on the electronic 
structure of Refi^), while applications can be found in 



edtp = —V • j, we obtain the expression for the charge 
current density j a (r,t) 



-Hli-T *)V Q * CT (r, t) + — A a (r, t)p(v, t) 

I m £ — / UlC 



(3) 

where TL denotes the hermitian part, and we assumed the 
vector potential to be directed along the a-direction. In 
the Coulomb gauge, the second term in Eq. (J3j) is the 
diamagnctic current, which we will drop hereafter, since 
its contribution to the conductivity is purely imaginary. 
The first term is called the paramagnetic current. For 
a discussion on their physical interpretation see e.gp^. 
The current expectation value is, within the Kubo linear 
response formalism, linked to the current-current corre- 
lation function. For the optical conductivity, we then 
find 



(4) 



where x(w) is the longwavelength-limit (q = 0) of 



B. The optical conductivity 

The optical conductivity tensor er Q/3 (q, ui) is defined as 
the linear response that relates the total electric field in 
the solid to the charge current density^ : 



X Q/3 (q,r) 



(j a (r,t))=Y,^(^^E f3 (r,t) 





(1) 



here a, (3 denote cartesian coordinates and (•) indicates 
the quantum mechanical expectation value. 

In the following we will derive an expression for the 
long wavelength limit (q = 0) of the real part of the 
conductivity tensor, which we shall refer to as the op- 
tical conductivity $i<j a P(uj). The starting point of the 
derivation is the fundamental Hamiltonian of the sys- 
tem, H = Hq + Hint, with Hq being the one-particle 
part, with the coupling to the (classical) light field via 
its vector potential A(r, t) : 



E 



/ d 3 r *t( r>t ) 



-l(^V + ^A(r,i)) 2 + y(r) 



(2) 

where we have chosen the Coulomb gauge, V • A = 0. 
V(r) is any one-particle potential. In practice, it will e.g. 
be the effective Kohn-Sham potential of density func- 
tional theory within the LDA. We emphasize that this 
notation in terms of the field operators, 'J, is still basis 
free. Hint contains, in our case, electron-electron interac- 
tions, involving only two-body terms. Then, it commutes 
with the charge density and is trivially gauge invariant^ 
Gauge invariance of the full Hamiltonian leads to 
charge conservation and, via the continuity equation 



-— <Tf(-q,T),/(q,0)) 



(5) 



which we have written in the imaginary time Matsub- 
ara formalism. So far, all quantities live in the spatial 
continuum. As announced, we shall make the connec- 
tion with many-body techniques that work in localized 
basis sets. At this point, however, we shall first de- 
velop the field operators in a Bloch-like basis, \I/(r, t) = 

EkL CT ( kL l r ) c kL CT ( r )- Hcrc > L = denotes or- 

bital (n,l,m) of atom 7 within the unit cell. c kL(T , c kLcr 
are the usual (discrete) creation and annihilation op- 
erators. The momentum sum runs over the first Bril- 
louin zone. Later on, we will switch to the Wannier-likc, 
real space basis |RX) = CX P (*kR) (ki) , to which the 
aforementioned notion of localization will apply. Here, R 
labels the unit-cell, which only in case of a one-atomic ba- 
sis is equivalent to the atomic position. This distinction 
will prove important later on. Taking the limit of long 
wavelengths, which in this context is the familiar dipole 
approximation, we find for the paramagnetic current 



j«(q = 0,r)=e 



k,LL'(T 



7, L V 
1>V ™ t-l. 



(rK La (r) (6) 



with the so-called Fermi velocity, or dipole matrix ele- 
ment 



"k.'a = -n(kL'\v a \kL) 



(7) 



where V a is the a-component of the momentum oper- 
ator. In the evaluation of q = correlation functions 
in infinite dimensions, thus in a dynamical mean-field 
spirit, vertex corrections are absent in the one orbital 
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case^iSi. In other words, electron-hole interactions effec- 
tively vanish, and the two-particle quantity, Eq. © , can 
be decoupled into the product of two one-particle Green's 
functions. In the following, we shall neglect vertex cor- 
rections, though the above statement is not valid for the 
cluster or multi orbital case. Indeed, one can show that 
only the elements of the Fermi velocity that are diagonal 
in orbital space have the odd parity with respect to mo- 
mentum that is required for the vanishing of vertex cor- 
rection in d = oo. Instead of advancing towards a more 
stringent two-particle formulation in the model case, it 
is our objective to strive after a formalism that accounts 
for the complexity of realistic state-of-the-art electronic 
structure techniques, such as LDA+DMFT. With this 
simplification, the derivation, after continuing to real fre- 
quencies {iLD n —> lu + i0 + ) , yields : 

v r J 

x trj^4(k,a/ +w)v a (k)A(k,w>/j(k)}( 8 ) 

Where tr stands for the trace over orbitals and A(k, ui) = 
— 1/71-3 [u + \i — H (k) — £(w)]~ and v a (k) are the or- 
bital matrices of the momentum-resolved spectral func- 
tions and the Fermi velocities, respectively. Ho denotes 
the one-particle Hamiltonian in the absence of the exter- 
nal field. Hence, in total, the conductivity acquires the 
well-known form of a frequency convolution of momen- 
tum resolved spectral functions, with the Fermi velocities 
modulating the amplitudes of the spectral weight. Due 
to the fact, mentioned earlier, that the interaction part of 
the Hamiltonian commutes with the charge density, the 
Fermi velocities are those of the non-interacting prob- 
lem Ho. The many-body physics only enters the spectral 
functions via the evaluation of the expectation value in 
the correlation function of Eq. (J5J). 

III. FERMI VELOCITIES 

The Fermi velocity matrix element, Eq. iJTj). is read- 
ily evaluated when e.g. working within a plane wave ba- 
sis set. Yet, many-body techniques, such as DMFT and 
its realistic extensions, that are geared at improving on 
local interactions in the spirit of the Hubbard model, 
necessitate the use of localized orbital sets, e.g. muffin 
tin derived orbitals, £/iVMT0 71 i 72 , or other Wannier 
functions 7 -^. While a computation of the full matrix el- 
ement, Eq. |(7J), is in principle still possible within these 
basis sets, it becomes rather tedious from the practical 
point of view. Therefore we have devised a handy ap- 
proximation, explicitly geared at the use with localized 
orbitals, which allows for a reliable calculation of optical 
properties at a rather low computational cos t 66 i . 

A. Lattice formulation : 



The Peierls substitution and its generalization 

In the above, we coupled the light field to the electronic 
degrees of freedom of the solid via the standard minimum 
coupling, and developed the continuous field operators 
into a basis, which led to the given Fermi velocities of 
Eq. ([7]). A different approach is to instead develop first 
the Hamiltonian in this basis, and to couple the vector 
potential directly to the site, or lattice operators in a 
way that verifies gauge invariance. Consider the Hubbard 
model 

H = E ^ L ^c ]La + n mt (9) 

ij,LL' ,a 

Then, the philosophy of the "Peierls substitution" 74 i 75 
approach is to add the following phase-factors to the lat- 
tice operators^ : c \ La -> c\ La exp (1% J RtL dr A(r, t)) . 

Here, R^l denotes the atomic positions. These can be 
separated into 

R, L = R, + Pl (10) 

where the former indexes the unit cell i and the latter 
the atom 7 within the cell, in the case of a multi-atomic 
unit cell. This seemingly trivial statement will lead to 
important terms in the Fermi velocities, which to our 
knowledge have so far not been considered. Equally, in a 
multi- atomic environment, the lattice position operator 
TZ can be defined as 

H = Y,KiLc\ L<r C iLc (11) 
iha 

When now supposing the interactions in Eq. ([9|) to be 
only of density-density type in the lattice operators (cf. 
footnote above), then the above phases only appear in 
the kinetic part of the Hamiltonian. When additionally 
assuming a slowly varying vector potential such as to 
approximate the integral in the exponent, the Peierls 
approach can also be seen as a substitution for the 
hopping amplitudes in the above lattice Hamiltonian : 
tf- L -> tfj L exp (i-t c A(t) (Rtf, - R jL )). We further re- 
mark that evidently the vector potential only couples 
to non-local hopping elements, i.e. within this approach 
intra-atomic transitions, = (n,l, 771,7)) ~~ * {hL — 

(71', I', m' , 7)) in the above notation, are absent. From 
the thus defined substitution, we can compute the cur- 
rent either by means of the continuity equation, or by a 
functional derivative of the Hamiltonian with respect to 
the vector potential. One finds 

j Q = e £ vt' a L ct L , a c kLa (12) 

LL',k,o- 

with the velocity 

vt' a L = ^tt; L (B? L ,-R? L )e-^-^ (13) 
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Using the separation Eq. (fit)]) , and introducing the usual 
Hamiltonian element H£ ' l = - tf- L e~ lk(R '- R i\ we 
find a generalized Pcicrls expression : 

vfa = I (dk a ^ L i(pt, p a L )nZ'A (14) 

The first term is the familiar Fermi velocity, given by the 
momentum derivative of the Hamiltonian. It contains 
hopping processes that take place between different unit 
cells While absent in the one- atomic case, 7 = 7', 
the second term, which to our knowledge is new, be- 
comes crucial, once calculations of realistic materials are 
performed. It accounts for hopping amplitudes between 
different atoms 7', 7 within the same unit cell. Indeed 
when considering e.g. a simple cubic one-atomic system, 
first in its primitive unit-cell, and then in a non-primitive 
unit cell that is doubled in the direction along which the 
momentum derivative is taken. Then it is the second 
term in Eq. (fT4)) that assures that the optical conduc- 
tivities of the two equivalent descriptions are the same. 
When using the derivative term only, the optical conduc- 
tivity comes out wrong, even qualitatively. 

The above formula is very handy, since the only in- 
tervening matrix clement is that of the Hamiltonian and 
its momentum derivative. No other integrals involving 
the LDA wave functions, which are cumbersome to han- 
dle in the chosen basis, occur in this case. We note that 
the above expression is hcrmitian. Yet, in general, it has 
no well defined parity with respect to the momentum k, 
even when assuming inversion symmetry of the Hamilto- 
nian. Only the elements that are diagonal in the atomic 
7-indices have the required odd parity that leads to the 
cancellation of vertex corrections in the limit of infinite 
coordination. 



B. Continuum formulation : 
Assessing the Peierls substitution 

In the preceding section, an expression for the Fermi 
velocity was deduced from the lattice formulation of the 
solid. This has to be contrasted to the proper matrix 
element, Eq. ([7]), that originates from the continuum de- 
scription. Now, a valid question is whether, and un- 
der which circumstances, the generalized Peierls velocity 
can be employed as a reliable approximation to the true 
dipole matrix element. In the appendix we show that the 
latter can actually quite naturally be split into the gen- 
eralized Peierls expression and a correction term that re- 
covers the full matrix clement. The impact of this supple- 
mentary term decreases with an increasing localization of 
the basis functions. Indeed, in the limit of strongly local- 
ized orbitals, the only missing terms are atomic transi- 
tions, which, as discussed above, are absent by construc- 
tion in the (generalized) Peierls approach, and have to be 
accounted for separately^. Therewith, the above derived 
expression of the Fermi velocity is particularly suited for 
use in Wannier function setups for compounds with d 



or f orbitals, which satisfy the requirement of localized 
orbitals. 

In the appendix, we give explicit expressions for going 
beyond the generalized Peierls approach, e.g. we derive 
a formula for including intra-atomic transitions within 
the current setup of localized orbitals. In the practical 
calculations for VO2 within a localized basis, however, 
we found these terms to be negligible as evidenced by 
the good agreement between the Peierls treatment and 
experimental findings. 



C. Downfolding of Fermi velocity matrix elements 
Upfolding of the downfolded response 

Many-body calculation for realistic systems often work 
in a downfolded setup. In other words, after a band- 
structure has been obtained from e.g. an LDA computa- 
tion, orbitals that are supposed to be subject to only mi- 
nor correlation effects are integrated out and linearized. 
These are typically high energy excitations, and thus 
the downfolding procedure is used to construct an ef- 
fective low-energy problem, which is simpler to be tack- 
led with a many-body approach. The linearization step 
preserves the Hamiltonian form of the one-particle part 
of the problem. Thereby the influence of correlation ef- 
fects beyond the one-particle band-structure of these or- 
bitals, and also the possible feedback on the others, are 
entirely neglected, and the high energy excitation spec- 
trum fixed to the Kohn-Sham eigenvalues within LDA. 
The many-body calculation thus lives in an orbital sub- 
space only, and all other orbital degrees of freedom re- 
main unaffected. The introduction of a double-counting 
term, which corrects for correlation effects already taken 
into account by the LDA, adjusts the center of gravity 
of the many-body spectrum with respect to the higher 
energy parts. See however e.g. Ref.— for a way how to 
include uncorrelated orbitals within the LDA+DMFT cy- 
cle. For a recent scheme to incorporate also self-energy 
effects of higher energy orbitals into an effective model 
see Refill. 

Although in the computation of the Fermi velocities, 
Eq. ([7]), only the "non-interacting" Hamiltonian enters, 
several complications occur, when it comes to deducing 
of optical properties from downfolded many-body calcu- 
lations : Not only are transitions from and to high energy 
orbitals truncated, but also the optical transitions within 
the block of low-energy orbitals acquire wrong ampli- 
tudes. This owes to the fact that, evidently, the compu- 
tation of transition matrix elements and the downfolding 
procedure do not commute. Here one has to distinguish 
between the effect on the full matrix element from that 
on the approximation of the Peierls velocity. Indeed, the 
orbitals of the downfolded system are in general less lo- 
calized than the ones of the original problem, and the 
Pcicrls approximation therewith is less accurate. For an 
instructive discussion on this subject see also££. 
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Here wc explain a simple strategy^ for the computa- 
tion of the optical conductivity, applicable to many-body 
electronic structure calculations that were performed us- 
ing a downfolded one-particle Hamiltonian. In line with 
the above remarks, this procedure yields better results 
than when computing the Fermi velocities directly from 
the downfolded system. The procedure is, moreover, not 
limited to the use of the Peierls approach. 

The central quantity to look at in this respect is the 
orbital trace of the matrix product of Fermi velocities 
and momentum-resolved spectral functions in Eq. ([5]) : 

tT^VkAk(uj')vkAk(u}' + w)j (15) 

Since the trace is invariant under unitary transforma- 
tions, the above can be written as 

tr(ulv k U k A k (Lj')ulv k U k A k (u}' (16) 

for arbitrary unitary matrices U k . In the case of a 
band-structure calculation (i.e. a vanishing self-energy, 
E = 0), we can chose these matrices such that they per- 
form the desired downfolding, i.e. both, the spectral func- 
tions A k = U^A k U k and the transformed Hamiltonian 
will acquire a block-diagonal form. In the following we 
shall distinguish between the low energy block "L" and 
the high energy block "H" . An LDA+DMFT calculation 
will add local Coulomb interactions only to the former, 
which will result in a self-energy that lives in the "L" 
sub-block, while the orbitals of the "H" block will remain 
unchanged from the many-body (DMFT) calculation. In 
other words, since both sub systems are disconnected, 
the block-diagonality of the spectra is retained through- 
out the calculation. 

The idea is now to compute the dipole matrix elements 
from the initial full system, and then to apply the same 
basis transformation that blockdiagonalizes the Hamilto- 
nian also to the velocities. 

Clearly the downfolding procedure is not exact, since 
it linearizes the impact of the high energy orbitals. When 
solving the system with the full, non-downfolded, Hamil- 
tonian, the matrices that block-diagonalize the full sys- 
tem would not be the same. They would even depend on 
the frequency us due to the dynamical nature of the self- 
energy. Yet, when granting the approximative validity of 
the downfolding as such, and assuming the U k to remain 
unchanged with respect to the initial band-structure, we 
can proceed further, and by specifying 

rrt tt f V x W\ 7 , „ / L \ 
v k = Kv k U k =i^ w f y 2 J, A k (us ) = Q H j, 

and A k (us' + us) = ( J (17) 

where the spectra of the L sector are taken from the 
many-body calculation. Then the above trace becomes 

LVxLVi + LWHW^ + HV 2 HV 2 + HW^LW (18) 



For transitions within the block of only correlated or- 
bitals, L , intervenes the Fermi velocity matrix V\, which 
is evaluated as the low-energy block of the unitary trans- 
formed matrix element of the full, i.e. non-downfolded 
system. The resulting velocity v k is thus different from 
the matrix element that is computed from the downfolded 
system. When using in particular the Peierls expression 
the momentum derivative of the unitary matrices U k lead 
to additional terms in the latter case. 

Moreover, with the above, a restriction to the low- 
energy block is not imperative. We can indeed calcu- 
late the complete optical response, including transitions 
from, to and within the high energy block^. The lat- 
ter is entirely determined by the band-structure calcu- 
lation. When comparing to experiments, this allows to 
assess whether it is only the relative position with respect 
to the low energy orbitals that needs an "adjustment" , 
which is connected to the double counting term, correct- 
ing for correlation effects already taken into account by 
the LDA, or whether correlation effects modify substan- 
tially the overall spectrum of downfolded orbitals. The 
latter can be brought about e.g. by non-negligible life- 
times, or shifts that depend on the individual orbital. 

We will refer to the above described scheme as "up- 
folding", since the downfolded orbitals are reintroduced 
for the sake of accounting for optical transitions from, 
into, and between them. 



IV. OPTICAL CONDUCTIVITY OF 
VANADIUM DIOXIDE - AN APPLICATION OF 
THE FORMALISM 

In a recent wor k 34 ' 35 , we used an analytical continu- 
ation procedure to calculate real-frequency self-energies 
from LDA+Cluster-DMFT d&tiM. This allowed for a 
better understanding of the impact of correlation effects, 
especially for the insulating phase of VO2. Indeed, it 
was revealed that correlation effects enhance the Peierls 
effect, while non-local fluctuations preserve the coherence 
of one-particle excitations. Moreover we successfully con- 
structed a static, yet orbital-dependent potential that, 
when added to the LDA Hamiltonian, reproduced the 
LDA+CDMFT excitation spectrum when artificially ne- 
glecting all lifetime effects. In our picture, Ml VO2 is 
not a genuine Mott insulator, and we referred to it as a 
"many body Peierls" insulator. 

Here, wc compute the optical conductivities of both, 
the metallic and the insulating phase. A comparison 
with experimental results will allow to further confirm 
the underlying LDA+CDMFT electronic structure cal- 
culation^, and its interpretatio n 34 ' 35 . Furthermore, this 
will also enable us to analyze the experimentally mea- 
sured intensities with a solid theoretical background. 

The LDA+DMFT optical conductivities are presented 
in the Figs. [1] and [3] for the metal and the insulator, re- 
spectively. In order to compute the optical conductivity 
also for high energies, we employ the upfolding scheme 
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detailed above. In the many-body Cluster-DMFT cal- 
culation^ all orbitals other than the vanadium t2 g were 
downfoldcd. The latter thus constitute the low energy 
sector, L, according to Eq. (jTTJ) . For the calculation of 
the Fermi velocities we use a larger Hamiltonian that 
comprises for the high energy part, H, in particular the 
vanadium e^ and the oxygen 2p orbitals, and, moreover, 
the oxygen 2s^. We sketchily write s,p,e^ in the graph- 
ics. When indicating that transitions are from s,p,eg into 
the t2 g orbitals, this mainly accounts for transitions from 
the occupied 02p into empty t2 S orbitals, since, e.g., the 
Cg to t-2g transitions are derived only from the little oc- 
cupied weight of e° character that stems from hybridiza- 
tions with occupied orbitals. 

When referring to the orientation of the electric field, 
or the light polarization, we use the simple monoclinic 
lattice as reference^. Since for the Peierls Fermi veloc- 
ity, Eq. (|A3j) , we perform the numerical derivative of the 
Hamiltonian on a discrete momentum mesh, not all direc- 
tions are accessible in a straight forward manner. Yet, 
the important polarizations, E \\ [001] and E _L [001], 
are capturable. In an experiment, the polarization is 
varied by choosing different orientations of the sample, 
or different substrates, which, in the case of thin films, 
favor different growth directions. Herewith, all orienta- 
tions that lie within the plane of the surface are probed, 
when using unpolarized light. In our calculations, how- 
ever, we evaluate the response of a single given polariza- 
tion only, without averaging over an ensemble of in-plane 
directions. 

As a comparison to our theoretical curves, we include 
results from three experiments that we already men- 
tioned in the beginning. We will display measurements 
on single crystals by Verleur et al£^, performed for differ- 
ent orientations of the sample. Moreover, recently, exper- 
iments were carried out on different types of thin films. 
The work of Okazaki et alX used thin films (T c w 290 K) 
with [001] orientation, i.e. for the electric field E _L [001]. 
Qazilbash et aZ.— on the other hand used polycrystallinc 
films with preferential [010] orientation (T c ~ 340 K). 
We now proceed with the presentation of our results for 
the individual phases. 



A. Rutile V0 2 - The metal 

In Fig. [1] we show, along with the mentioned experi- 
mental data, the theoretical optical conductivity of ru- 
tile VO2, which we obtain for the different light polar- 
izations as indicated. As one can see, already the three 
experiments yield quite distinguishable spectra. The dif- 
ferences may point to a polarization dependence, but one 
cannot rule out an influence of the sample type and the 
means by which multiple reflections at the sample sub- 
strate were treated in case of the thin films. Indeed, in the 
case of rutile VO2, x-ray experiments^ witness a rather 
isotropic response. The different measurements on single 
crystals^ also evidence a quite uniform conductivity up 
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FIG. 1: (color online) LDA+CDMFT optical conductivity 
of the rutile phase of VO2 for the indicated polarizations 
([aab] = [0.85 0.85 0.53]). The velocity matrix elements were 
calculated using the scheme of Section MI CI Beyond the 
t2 9 orbitals this calculation includes in particular the Ve^ 
and 02p orbitals. Experimental curves from (a)— (single 
crystals, orientation as indicated), (b)-^ (polycrystalline film 
(T c « 340 K), preferential orientation E _L [010], T=360 K), 
and (c)i (thin film (T c « 290 K), E 1 [001], T=300 K). 




FIG. 2: (color online) LDA+CDMFT optical conductivity of 
rutile VO2 for the [001] polarization. Shown are the differ- 
ent orbital transitions according to their energy sector, see 
Eq. (|18p . The contributions are additive and sum up to the 
total conductivity. For details see Section l"!!! CI Experimental 
curves, as above fromi^. 



to 4 eV. The polarization dependence of the theoretical 
conductivity is found to be rather small, too, which is also 
in line with our previous statement^ that the t2 5 self- 
energy shows no particular orbital dependence. Thus, in 
theory, the metallic Drudc-like response is made up from 
ai s and Cg density near the Fermi surface. 

At higher energies, beyond the Drude-like tail, further 
inter- "band" intra-t2 S transitions occur. Yet, the optical 
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response is rather structureless up to 2 cV. At this en- 
ergy, however, we already expect the onset of oxygen 2p 
derived transitions. In order to elucidate the origin of the 
spectral weight of this region in greater detail, we plot in 
Fig- ID the optical conductivity resolved into the different 
energy sectors, according to Eq. (fT8|) . Since the 02p and 
the orbitals were part of the downfoldcd high energy 
sector, their position, within our scheme, is frozen to the 
LDA result (see e.g. the band-structure in22). Therefore 
transitions from the 02p orbitals into the t 2£( ones start, 
as expected, at around 2 eV. We remark that the po- 
larization dependence for the oxygen derived transitions 
agrees very well with the single crystal experiments 4 ^ up 
to 4.5 eV. Transitions from the t2 g orbitals into the e^ set 
in later, at around 2.5 eV, and are rather small in magni- 
tude. The 02p to transitions appear at the expected 
energies, but they are too low to be seen in Fig. [2] 

Overall, the LDA eigenvalues seem to give a rather 
good description of the e a g and 02p orbitals, since the 
agreement with experiment is reasonably accurate, as 
was qualitatively noticed already in previous LDA op- 
tics calculations^. When looking at photoemission rc- 
sult a 21 i 46 , one remarks that the on-set of the oxygen 2p 
is compatible with the LDA, yet, their center of gravity 
is shifted to slightly higher binding energies in the ex- 
periment. As to the e^ orbitals, it is conceivable, when 
resorting to x-ray experiment a 21 i 46 as a reference, that 
they appear at a little larger energies and with a smaller 
bandwidth than within the LDA. Of course both com- 
parisons are somewhat indirect, due to the occurrence 
of matrix elements and other effects in the experiments. 
Yet, we emphasize that the rather incoherent nature of 
the t2g weight in the spectral function^ is far beyond any 
band-structure technique, which is why the optical con- 
ductivity in the 2.5 to 4.0 eV region, derived from 02p 
to t'2 g transitions, comes out too large in LDA^ when 
comparing to the experiment of Ref.— , while we find a 
good agreement for the LDA+CDMFT conductivity. 

At this point, we can only speculate on the origin 
of the shoulder and peak structure seen in one of the 
experiments 4 ^ at 2.5 eV, and 3.0 eV, see Fig. [51 It 
seems conceivable that it stems from %2 g to 02p transi- 
tions, rather than from e^ contributions. Attributing the 
humps to distinct 02p to a l9 or e^ transitions is cum- 
bersome, mostly due to the structure of the numerous 
oxygen bands. When looking at the momentum-resolved 
optical conductivity (not shown), one realizes that 02p 
to Cg transitions start for most of the k-regions at lower 
energies than transitions into the ai g . 



B. Monoclinic VO2 — The insulator 

We now come to the optical spectra of the monoclinic 
phase. In Fig. [3]we show our theoretical LDA+CDMFT 
results, again, in conjunction with the three experi- 
ments ^ 44 ^ 48 . As was the case for the metallic phase, 
the latter yield varying results. While the optical gap 
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FIG. 3: (color online) LDA+CDMFT optical conductiv- 
ity of the Ml phase of VO2 for the indicated polarizations 
([aab] = [0.84 0.84 0.54]). The velocity matrix elements were 
calculated using the scheme of Section IIII CI Experimental 
curves from (a)^ (single crystals, orientation as indicated), 
(b)i (thin film, E _L [001], T=280 K), and (c) 4 ^ (polycrys- 
talline film, preferential orientation E _L [010], T=295 K). 
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FIG. 4: (color online) Optical conductivity of the Ml phase 
of VO2 for the indicated polarizations ([aab] = [0.84 0.84 0.54]) 
when using the effective band-structure^ "LDA+A" . Exper- 
imental curves from (a)^i (single crystals, orientation as indi- 
cated), (b)i (thin film, E 1 [001], T=280 K), and (c) 4 ^ (poly- 
crystalline film, preferential orientation E _L [010], T=295 K) 



is roughly 0.5 eV in all cases, the higher energy response 
is markedly different. Not only the amplitudes, but also 
the peak positions differ considerably. Yet, as a matter 
of fact, in the current case of Ml VO2, a sizable polariza- 
tion dependence is expected from the structural consider- 
ations mentioned above. Indeed our calculation suggests 
a noticeable anisotropy in the optical response, which is 
congruent with the experimental findings. 

Before interpreting the results, however, we find it 
instructive to also compute the conductivities using a 
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somewhat simpler approach : In Ref i ' 35 we deduced 
from the LDA+CDMFT self-energies an effective static, 
yet orbital-dependent, one-particle potential that repro- 
duced the many-body excitation spectrum, which arises 
when neglecting all life-time effects. Therewith all corre- 
lation induced energy shifts are captured, whereas the co- 
herence of the excitations remains infinite. This is equiv- 
alent to the use of a scissors operator, albeit one that does 
not simply widen the gapS 3 -, but selectively shifts the 
one-particle excitations that mediate the dimerization. 
We thus replace in a full orbital LDA Hamiltonian the 
Kohn-Sham eigenvalues of the t2 g orbitals by the ones to 
which the additional potential was applied, and label the 
results "LDA+A". For further details see 3 ^. There- 
with, we do not have to invoke the upfolding scheme of 
Section [mCl 

What this theoretical conductivity is missing are the 
life-time effects encoded in the imaginary part of the 
LDA+CDMFT self-energy. These were found to be 
small, yet not entirely negligible 3 - 5 -. Fig. [3] displays our 
result, again along with the experimental curves. 

When looking first at the optical conductivity that re- 
sults from the effective band-structure 3 ^ 5 -, "LDA+A" , 
Fig. 21 we find all experimental polarization tenden- 
cies reproduced : Consistent with Verleur et ah— , the 
E || [001] conductivity is lower than the E L [001] one 
at energies up to 1.5 eV, after which the c-axis response 
develops a little maximum of spectral weight in both, ex- 
periment and theory. At energies of 2.35 eV— or 3.0 eV— 
the experimental conductivity with E |j [001] compo- 
nents evidences a narrow peak. In the calculation this 
is prominently seen at 2.75 eV. When looking at our ef- 
fective band-structur o 34 ' 35 , it seems plausible that these 
transitions stem from ai 9 bonding to anti-bonding or- 
bitals. The peak is indeed very narrow for an inter-band 
transition, but in our picture this is simply owing to the 
fact that the B,\ g anti-bonding excitation does exhibit an 
almost dispersionless behavior— i 3 - 5 -. However, already in 
this frequency region we expect transitions that involve 
the oxygen 2p orbitals, as will be detailed below for the 
LDA+CDMFT conductivity. 

At still higher energies, the E || [001] response is again 
lower than for the perpendicular direction in both, ex- 
periment and theory. The overall congruity with exper- 
iments further corroborates the validity of our effective 
band-structure picture for spectral properties and there- 
with strengthens our interpretation of the nature of the 
insulating phase of VO2 as a realization of a "many-body 
Peierls" stated 5 -. 

Now we compare this simplified approach with the full 
LDA+CDMFT conductivities of Fig. H We instantly 
realize that the LDA+CDMFT response for the t2 S or- 
bitals is damped and therewith less structured, which 
was clearly expected. The small underestimation of the 
optical gap is probably owing to the elevated tempera- 
ture at which the LDA+CDMFT quantum Monte Carlo 
calculation was performed^ 3 -. 

To shed further light on the structure of the re- 
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FIG. 5: (color online) LDA+CDMFT optical conductivity of 
Ml VO2 for the [001] polarization. Shown are the different or- 
bital transitions according to their energy sector, see Eq. (|18l) . 
The contributions are additive and sum up to the total con- 
ductivity. For details see Section [HI CI Experimental curve, 
as above fromi^. 



sponse, we resolve in Fig.[5]the contributions to the [001] 
LDA+CDMFT conductivity into their respective energy 
sectors, according to Section IIII CI From this we first 
infer that the slight upturn, seen for this polarization 
beyond 1.5 eV in the LDA+CDMFT conductivity is in- 
deed derived from transitions within the t2 g manifold, for 
oxygen contributions only set in at around 2.0 eV. 

Besides, the prominent peak in both, the experimen- 
tal and LDA+A conductivity with E || [001] polarization 
that we attributed above to ai g -ai s transitions, is largely 
suppressed and only faintly discernible as a weak shoul- 
der, when comparing with the other polarizations. 

As an explanation for this difference between experi- 
ment and the approach of the one-particle potential A on 
the one hand, and the LDA+CDMFT result on the other, 
we forward the occurrence of sizable life-time effects in 
the LDA+CDMFT electronic structure calculation. In- 
deed the &i g spectral weight in the corresponding t2 9 
LDA+CDMFT spectral function is not sharply defined 
and extends over more than 2 eV, and is only barely dis- 
cernible in the total, orbitally traced, spectrum 3 - 3 -. When 
thinking of the conductivity in simple terms of density- 
density transitions, it is perfectly conceivable that the 
&i g -&i g response eventuates only in a tail of spectral 
weight (as seen in the energy sector resolved conductivity 
in Fig. [5]) and not in a well defined peak. Having said 
this, and referring to the experiments, one might thus 
conjecture that these life-time effects are still overesti- 
mated in the LDA+CDMFT calculation. Moreover, we 
stress again that the many-body electronic structure was 
computed at high temperatures^ 3 -, which will lead to an 
overestimation of the temperature induced part of the 
broadening^ 3 -. 

Finally, we remark that despite all differences in the 
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experimental data, they reveal (maybe apart from the 
single crystal for [001] polarization) a common global 
tendency namely that, when going from the metal to 
the insulator, low frequency spectral weight is transfered 
to higher energies. Indeed, for a given polarization, the 
Drude-like weight that the insulator is lacking at low en- 
ergies must be recovered, as requires the f-sum ruleZ^. 
This condition is met at 5.5 cV in one experiment^, 
while in the other— an overcompensation appears already 
at energies beyond 3.5 eV. Theoretically, when using the 
LDA+CDMFT conductivity, we find values of 3.73 eV 
and 4.35 eV, for the [110] and [001] direction, respec- 
tively. 



C. Comparison with simpler approaches 

In this section we shall briefly show that all the trouble 
with the Fermi velocity is worth the effort. Therefore we 
plot in Fig. [6] for the case of Ml VO2 a comparison of our 
full scheme, which proved to yield quantitative results, 
with two simplified calculations. These differ from the 
full scheme only in the way how the Fermi velocities, i.e. 
the transition amplitudes are treated. We restrict the 
discussion to the t2 g response. 

To illustrate the effect of the downfolding of orbitals on 
the matrix elements, we have computed the optical con- 
ductivity when applying the generalized Peierls formula 
on the downfolded Hamiltonian. As we can see, the re- 
sulting curves differ considerably from those using the 
upfolding scheme. In particular, the absolute value for 
some polarizations is way off with respect to experiment. 

It has become a popular approximation to entirely ne- 
glect Fermi velocities in the computation of optical prop- 
erties. Therewith the conductivity is a simple convo- 
lution of momentum-resolved spectral functions. As a 
consequence inter-band transitions are omitted, since the 
Fermi-vclocities are simple unit matrices. Especially in 
the realistic context this is a severe oversimplification. 

Moreover, by construction, there cannot be any orbital 
dependence in the conductivity, while, as evidenced from 
the experiments, this is clearly an important issue for Ml 
VO2. Also, the magnitude of intra-band transitions is not 
properly accounted for. In fact, the absolute value of the 
response is not well defined. In order for this approach 
to yield a comparable magnitude, we arbitrarily choose 
a pre-factor : v^ L = 2r ^LL', with r being the Bohr 
radius. As can be inferred from Fig. [6l the resulting 
peak structure of the optical conductivity is wrong. 

Finally, we also compute the optical conductivity when 
neglecting the multi-atomic correction term in the veloc- 
ities, i.e. using only the derivative of the Hamiltonian. 
As was the case for the velocities of the downfolded case, 
only for one polarization does this yield a reasonable re- 
sult. 
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FIG. 6: (color online) LDA+CDMFT t 2g optical conductivity 
of the Ml phase for the indicated polarizations ([aab] = [0.85 
0.85 0.53]) and different Fermi velocities. The "upfolded" 
curves correspond to our full scheme. The "downfolded" data 
computes the Fermi velocities from the downfolded Hamilto- 
nian. "no velocites" refers to a simple convolution of spectral 
functions without transition amplitudes. The "dH/dk" curves 
correspond to neglecting the multi-atomic generalization in 
the velocities. 



V. CONCLUSIONS 



In conclusion we presented a versatile scheme for the 
calculation of optical properties of correlated materials. 
Geared at the use with a localized basis set, we devised 
a realistic extension of the Peierls substitution approach. 
Moreover, we developed means to incorporate transitions 
that involve high energy orbitals that were downfolded in 
the many-body treatment of the electronic structure. 

As an application, we evaluated the optical conductiv- 
ity of VO2 for both, the metallic and the insulating phase. 
While the metal is characterized by a rather isotropic re- 
sponse, the insulator revealed a noticeable polarization 
dependence. The agreement with experiments is overall 
satisfying. The high energy conductivity is reasonably 
described when using the LDA band-structure for high 
lying orbitals. The LDA+CDMFT many-body calcula- 
tion for the t 2£( orbitals correctly describes the low-energy 
behavior. In the rutile phase it accounts for life-time ef- 
fects within the t2 ff orbitals and therewith also for the 
damping of oxygen to t2 9 transitions with respect to LDA 
results. In the insulator, it allowed for a genuine repro- 
duction of the experimental t2 S response, capturing in 
particular the polarization dependence over a wide en- 
ergy range. The congruity of experiment and theory for 
the t2 g spectral weight can be interpreted as corroborat- 
ing the validity of the underlying many-body calculation 
for the electronic structure along with its interpretation. 
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APPENDIX A: CONTINUUM FORMULATION 
OF THE TRANSITION MATRIX ELEMENTS 



Derivation of the generalized Peierls velocity in 
the continuum 



Starting from the general Fermi velocity, Eq. ([7]) , which 
originated from the continuum formulation, we here red- 
erive the generalized Peierls expression as an approxi- 
mation. The correctional terms contain all intra-atomic 
transitions, that were completely lacking in the lattice 
theory, as well as contributions that are owing to the 
spatial extensions of the wave functions in the solid. 

Using V a = — im/H [JZ a , H ], the element of Eq. ([7]) 
can be written 

R.R' 



where we have chosen to condense everything into two 
different terms. The first one obviously is 



d 3 r 



(R'L'|ft Q |r)(r|H |RX) 



<R'L'|Ho|r)(r|ft Q |RX) 



(Al) 



It is important to note, that here the position operator 
lZ a is defined in the continuum. Its effect in the posi- 
tion representation is (r|7vL a |RL) = r a x RL {r). This is to 
be contrasted to the discrete lattice version of Eq. (fTTjl . 
Moreover, one has to make a clear distinction between 
the continuous space variable r and the discrete unit cell 
label R. In the (unphysical) limit of completely localized 
wave functions, (r|RX) = Xrj^ 1 ") ~ ^( r — ^-l), this dis- 
tinction is relaxed, and we recover the expression Eq. fTJJ 
of the Peierls approach, as we shall sec. Like in the lattice 
case, we split the atomic positions into R" = R a + pi- 
Then the above becomes 

I(kL'|P a |kL> = -~ Ve ik ( R - R ')x 
m fi N 

R.R' 

(R' a -R a + pf, - p2)(R'£'|H |R£) + 
j d 3 r r a [(R'L'\r + R' + p v ) (r + R + p v |H |RL) 
-<R'L'|H |r + R + p L )(v + R + p L \RL)] }(A2) 



i(p2,-pg)(kL'|Ho|kL> (A3) 



KI: <ki ' |Ho|ki> 



which is exactly the generalized Peierls expression, 
Eq. HU). The merit of the Peierls approximation, in 
particular in realistic calculations, is its apparent sim- 
plicity. Indeed no matrix elements other than the Hamil- 
tonian need to be calculated. The latter is a quantity 
that is anyhow required for a many-body calculation, and 
one thus has only to perform the directional momentum 
derivative^ From the discussion of the Peierls substitu- 
tion above, it is clear that the second term in Eq. (|A2|) ac- 
counts on the one hand for all atomic transitions (R = R' 
and 7' = 7), yet it also contains contributions that arise 
from the fact that we started from a continuum formula- 
tion. In other words, the spatial extensions of the wave 
functions lead to inter-atomic, 7' 7^ 7, corrections, owing 
to their finite overlap. Yet, a direct evaluation of these 
terms is an intricate undertaking, since it involves the cal- 
culation of many integrals. Therefore, it is a valid ques- 
tion whether the generalization of the Peierls approach, 
as such, already gives a reasonable approximation, with- 
out considering the terms beyond it, and if so, under 
which circumstances. Though the regrouping of terms 
into the Peierls expression and the rest was guided from 
the lattice considerations, it might still seem somewhat 
arbitrary. The next section however reveals that the intu- 
ition of an increased validity of the Peierls approach with 
a better localization of the involved orbitals is actually 
warranted. 



2. The Peierls substitution as the localized limit 

In the following, we will make consecutive approxima- 
tions regarding the extension of the orbitals, which lead, 
step by step, to more simplified correction terms to the 
Peierls expression Eq. (|A3|) , which one might endeavor to 
take into account in an actual computation. Moreover, 
these approximative steps will rationalize the identifica- 
tion of the Peierls term as the leading contribution to the 
Fermi velocities in the considered setup. 

By assuming well localized orbitals, we thus proceed 
to cut down the expression in question to the predom- 
inant terms, which will be given by the integrals that 
involve wave functions that have a large overlap^. In- 
deed, we show that in the limit of perfect localization (or 
equivalently in the limit of large atomic separation) the 
only surviving transition elements are given by the intra- 
atomic contributions, that were missing in the Peierls 
formulation. Using (r + R|RL) = Xrx(i" + R) — Xol( t )i 
the terms beyond the Peierls ones can be put into the 
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form 



I 1 



£ e* k ( R - R '> / d 3 r r Q £ (A4) 

R.R' ^ A.L" 

XoL'(r + PL')XoL-(r + R' - A + ^)(A£"|Ho|RL) 
-x5L-(r + R - A + p L )xoL(r + Pl)<R'£'|H |AL"> 



In this formula the origins of all intervening wave func- 
tions lie within the same unit cell, labeled "0"—. In a 
first step, the assumed localization of the involved or- 
bitals makes it reasonable to identify important terms in 
the sum as those, where the arguments of the wave func- 
tions also lie within the same unit cell, i.e. A = R' in the 
first term and A = R in the second one. We note that 
within this approximation, only the Hamiltonian clement 
depends on the unit cell labels R and R , and we can thus 
directly perform the Fourier transformation, yielding 



T.'l 



(A5) 



Xol> (r + PL' )Xol» (r + PlO <kL" | H |kL) 
-XoL"(r + PL)x 0L (r + p L )(kL'\U \kL") 



This means that the entire momentum dependence, in 
this approximation, is carried by the Hamiltonian. The 
complexity of the occurring matrix elements of the po- 
sition operator has been considerably reduced. In the 
one-atomic case, i.e. 7 = 7' = 7", and when using the 
short-hand notation R^' = (0L\TZ a \0L') , Ho L '(k) = 
(ki|Ho|kL'} we simply have 



- [R«,o,Ho(k)] M 



(A6) 



This is reminiscent of the relation 1/mV = —i/h [1Z, Ho], 
which we used in the beginning. Here however inter- 
vene on-site matrix elements rather than the full position 
operator. Indeed these elements, R^ L , are well known 
in atomic physics : They give the usual amplitudes for 
atomic dipolar transitions : The angular part of the inte- 
gral will produce the corresponding dipole selection rules 
(Al = ±1, Am = 0, ±1) via Clebsch-Gordon coefficients 
(see e.g>^), when, as we have assumed, the wave func- 
tions have a well-defined angular momentum (7, m). Con- 
trary to the atomic case, however, the Hamiltonian is 
momentum dependent, owing to the fact that, though 
regarding atomic transitions, the "atom" is here embed- 
ded in a solid. Also, the above term reminds the form of 
the multi-atomic correction term in Eq. (|A3[) . only that 
there occurred fixed atomic positions p 7 , which commute 
with the Hamiltonian, which is why in Eq. (|A3[) only the 
non-local terms 7 ^ 7' appear. 

Coming back to the multi-atomic case, we have to 
make a further approximation in order to obtain an ex- 
pression containing atomic transitions only. Yet, the 



shifts in the wave functions of Eq. (|A5|) can be treated 
analogous to the unit cell coordinates : Indeed Xol ( r + 
Pl) is centered around the position of atom 7. When, for 
the sake of clarity, we rename Xoi( r ) = Xoi( r + Pl) we 
find 



/ d 3 rr Q £ 

J Til 



(A7) 



X6L'(r)XoL»(r + p L - - p L ")(kL"\H Q \kL) 



-XoL»(r + PL - PL")XoL(r)(kL'|H |kL"> 



From this expression it is plausible, that for localized 
orbitals atomic transitions (7" = 7' and 7" = 7, re- 
spectively) arc in fact predominant. When restraining 
ourselves to these cases, we thus drop entirely the cor- 
rections to hopping processes that stem from the finite 
extensions of the wave-functions and end up with 



- J d 3 r r a XoL<(r)XoL»W<kL"|H |kL) 

- 7 E^L"(r)XoL(r)(kL'|H |kL") 

"=7' 

(OL'\TZ a \OL"){kL"\H \kL) 

L" 

y 

(0L"|ft Q |0L)(kL'|H |kL") 



(A8) 



Here, only the terms in which the Hamiltonian element 
is diagonal in the atomic index 7 can be written in the 
form of a commutator, as was the case in the one-atomic 
case in Eq. (|A6|) . In total, the above term contains intra- 
atomic transitions only. These were completely missing 
in the Peierls approach, as explained above. Under the 
assumptions on the localization of the involved orbitals, 
the Peierls term, Eq. (|A3[) , thus turns out to be the most 
important contribution to the Fermi velocity. Intuitively, 
this approach is thus particularly suited for systems in 
which e.g. 3d or 4f orbitals play an important role, since 
these verify the request of a high degree of localization. 

Another point, worth noticing, is the fact, that, when 
using the Peierls approximation, the result of the conduc- 
tivity is actually basis dependent. It is the momentum 
derivative of the Hamiltonian, which constitutes the first 
term in Eq. (|A3|) . that transforms evidently differently 
than the Hamiltonian itself. Obviously this is an artifact 
of the approximations from which the Peierls expression 
eventuated. Since, however, the term corresponds to the 
limit of perfect localization, it is expected to still yield 
reasonable results for orbitals that are short-ranged. We 
will come back to this in the next paragraph, in the con- 
text of Fermi velocities for downfolded Hamiltonians. 
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Improvements to the above approximations are obvi- 
ous: One could e.g. take into account elements contain- 
ing nearest-neighbor wave functions within the same unit 
cell, or even account for wave functions centered in differ- 
ent unit cells. An evaluation of these terms in principle 
allows for a more quantitative assessment of the quality 
of the Peierls term. However, the matrix elements that 



one needs to evaluate are numerous and more complex 
since they explicitly involve various wave functions. We 
stress again, that these terms are inter-atomic corrections 
to the Peierls term, while the intra-atomic contributions 
are absent in the Peierls formalism by construction, and 
given by Eq. l[A"8]) . 
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